Electron-phonon vs. electron-impurity interactions with small electron bandwidths 
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It is common practice to try to understand electron interactions in metals by defining a hierarchy 
of energy scales. Very often, the Fermi energy is considered the largest, so much so that frequently 
bandwidths are approximated as infinite. The reasoning is that attention should properly be focused 
on energy levels near the Fermi level, and details of the bands well away from the Fermi level are 
unimportant. However, a finite bandwidth can play an important role for low frequency properties: 
following a number of recent papers, we examine electron-impurity and electron-phonon interactions 
in bands with finite widths. In particular, we examine the behaviour of the electron self energy, 
spectral function, density of states, and dispersion, when the phonon spectral function is treated 
realistically as a broad Lorentzian function. With this phonon spectrum, impurity scattering has a 
significant non-linear effect. 

PACS numbers: 71.10.Ay, 71.20.-b, 63.20.Kr, 72.10.Fk 
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I. INTRODUCTION 

In the last few years, a number of authors have investi- 
gated the effect of a finite Fermi energy on the properties 
of an interacting electron system 0, S S S IE 0- 
In the past, with the understanding that only properties 
near the Fermi level were important, the Fermi energy 
was assumed to be infinite, to facilitate integrations over 
electronic energies 0. This approximation was believed 
to have negligible impact on properties concerning states 
near the Fermi level. Alexandrov et al. Q first noted that 
accounting for a finite bandwidth leads to a significant 
modification of the electron density of states (EDOS) 
near the Fermi level; they interpreted this as a "polaronic 
collapse" of the electron band. From our point of view 
the important observation is that the EDOS actually ex- 
periences a significant alteration with electron-phonon 
coupling. In this paper, following the work presented in 
Refs. U0, we wish to investigate further the impact of 
a finite density of states when both impurity scattering 
and electron-phonon coupling are present. Specifically 
we will examine the impact of a different bare EDOS (a 
Lorentzian model), and the impact of a more realistic 
broad phonon spectrum (previous work focussed on an 
Einstein spectrum). 

The paper is organized as follows: in the following 
section we briefiy review the formalism to calculate the 
electron self energy. The formulation we use allows for 
several properties to be obtained explicitly and analyti- 
cally. Section III deals briefly with only electron-impurity 
scattering. Significant changes occur at a very basic 
level because the bandwidth is finite. For example, the 
self energy acquires a real part 0. Section IV shows 
a number of results for combined electron-phonon and 
electron-impurity scattering, where both a delta-function 
and broadened Lorentzian form is used for the phonon 
spectrum. We also make a direct comparison with the 
results of Ref. 0| to show the model dependence of the 



analysis. We close with a Summary of our results. 



II. THE FORMALISM 

The self-consistent Migdal approximation has been ex- 
plained before in the cited references, so the reader is 
referred to them. For an arbitrary electron-phonon spec- 
tral function, a'^F{v), the resulting equations for the self 
energy, E(aj -I- i5), are as follows, at zero temperature: 
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Here, the electron Green function, G{e,Lu -\- iS), is given 
by the Dyson equation: 



G{e,uj + iS) = l/[uj + iS 



S(w + i(5)], 



(3) 



where the momentum dependence is summarized by the 
energy e. The parameter l/r characterizes the impurity 
scattering strength, and has units of a scattering rate. 
The bare electron density of states (EDOS) is given by 
No{e), while the renormalized density of states (RDOS) 
is given by N(lo). The former function in principle comes 
from a band structure calculation; we will simply model it 
either as a constant with finite bandwidth or a Lorentzian 
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[3|. The RDOS is a product of our calculation, and is 
given by 



iV(w) 



de No{e)A{e,Lo), 
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where A(e, 



— 1 



^ImG(e, Lu+iS) is the electron spectral 



function. For the two simple models for the bare EDOS, 
N{uj) can be determined analytically in terms of the elec- 
tron self energy 0. Thus, these equations constitute a 
set of self-consistent equations for the electron self en- 
ergy, from which all single electron properties can be cal- 
culated. As already mentioned, for the electron-impurity 
scattering, the process is characterized by a single num- 
ber, 1/t. On the other hand, the electron-phonon inter- 
action is specified by an entire function, a^Fiy), which 
we will sometimes take to be an Einstein spectrum, in 
which case two parameters are required, the frequency, 
and the "strength", as measured by, say, the mass en- 



hancement parameter, A = 2 Jp°° dv 



a^F(v) 



To specifi- 



cally examine the impact of broadened spectra, we will 
utilize the truncated Lorentzian shape defined in Ref. '3. 

Note that an alternative formulation first determines 
the relevant Green function and self energy on the imag- 
inary axis; then a set of self-consistent equations can be 
utilized to analytically continue the result to just above 
the real axis 0. This method was used in Ref. (00)) 
is particularly useful at finite temperature. 



III. ELECTRON-IMPURITY SCATTERING 

Eqs. H1I4|) require iterative self-consistent solutions. 
Analytic results are not possible, except in some limits. 
For example, with impurity scattering only, and with a 
bare EDOS given by a Lorentzian with half- width D/2 = 
ep, the self energy is given by 




(5) 



For an infinite bandwidth the impurity self energy is pure 
imaginary, '^{uj + i6) = finally an approximate re- 

sult is obtained by using the first term in Eqs. (|1I2|I with 
the bare Green function on the right hand side. These 
results are illustrated in Fig. 1. Similar results are ob- 
tained for a constant bare EDOS with finite bandwidth 
. For the imaginary part of the self energy (which Eq. 
El indicates is proportional to the density of states when 
only impurity scattering is present) the main effect of the 
finite bandwidth is to reduce the self energy to zero be- 
yond energies corresponding to the electronic states. The 
real part (Fig. la) experiences the biggest change, since it 
is zero at all frequencies if the limit of infinite bandwidth 
is taken first. Fig. lb illustrates that self-consistency im- 
pacts the low frequency behaviour in particular. Note 
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FIG. 1: Real (a) and imaginary (b) parts of the electron self- 
energy, with an electron-impurity scattering strength 1/r = 
1 meV. The results for different bandwidths are shown, for 
both the fully self consistent result (E) and for the non-self 
consistent result (Eq). The legend for both parts is given 
in (a). For infinite bandwidth the real part is zero, and the 
imaginary part is a constant at — l/(2r) (not shown). Note 
that full self consistency alters the scattering rate at the Fermi 
level. 



that while the intercept of the imaginary part of the self 
energy retains the physical interpretation of a scattering 
rate, the slope of the real part is no longer related to an 
effective mass. This will be pertinent when both electron- 
impurity and electron-phonon scattering are present, as 
already pointed out in Ref. 0. 

The renormalized EDOS is already shown (to within 
a minus sign and a constant factor — see Eq. Q) in 
Fig. lb. Note the renormalization at the Fermi level 
for the fully self-consistent results. From Eq. 01 we 
know that the renormalization of the EDOS is a convo- 
lution over the spectral function. Fig. 2 shows the spec- 
tral function at the Fermi level. As expected, the fully 
self-consistent calculation (Fig. 2a) reflects the renor- 



3 




> 

E, 

0) 
DC 




(a) 



(o(meV) 



60 80 
o)(meV) 



1/x(meV)=0 

5 

10 

20 

50 

_^ 75 -- - - 

100 120 140 



3 
< 



0.6 - 



0.4 



0.2 





D(meV)=5 




10 




20 




50 . 






(b) 





0.5 



1 

(o(meV) 



1.5 



FIG. 2: The full self consistent (a) and non-self consistent (b) 
electron spectral function for various values of the bandwidth, 
including the infinite bandwidth case. This latter result is of 
Lorentzian form, whereas all others deviate from this simple 
form. Note that the non-self consistent results (b) can be very 
non-Lorentzian for small bandwidths. 
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FIG. 3: Real (a) and imaginary (b) parts of the electron self- 
energy, for a moderate electron-phonon coupling strength, 
A = 1. The phonon spectrum is a delta function at fre- 
quency ve ~ 10 meV. Results are shown for one band- 
width, D — Wve , and for various electron-impurity scattering 
strengths, as shown. Note that the logarithmic singularity 
in (a) remains for all impurity scattering strengths. Over- 
all, the "sharp" features due to the Einstein oscillator remain 
throughout, independent of the impurity scattering rate. 



malization whereas the non-self-consistent one does not 
(Fig. 2b). The latter shows that the spectral shape can 
become very non-Lorentzian when the Fermi energy be- 
comes quite small. The former shows that it can be quite 
difficult to unravel the impurity scattering strength from 
the details of the electronic band structure (and here we 
are simply varying the bandwidth). As we shall see, the 
contributions from electron-phonon scattering generate 
even more ambiguity. 



IV. ELECTRON-PHONON PLUS 
ELECTRON-IMPURITY SCATTERING 

(A) EINSTEIN PHONON 

Results with just the electron-phonon interaction have 
already been presented 0, so we will focus on the two 
effects in combination. Figs. 3a and 3b show the real 
and imaginary parts of the self energy as the electron- 
impurity scattering strength is varied. For the electron- 
phonon scattering, we have used an Einstein spectrum 
with frequency i^e — 10 meV, and strength A = 1. Our 
results are similar to those of Ref. 4 except that our real 
part has much sharper features than theirs, presumably 
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FIG. 4: The dispersion corresponding to the parameters of 
Fig. 3. As mentioned with respect to Fig. 3, sharp features 
remain in spite of the impurity scattering. 



because ours is at T = (theirs is at a small non-zero 
temperature). It is also evident from our method of so- 
lution (see Eqs. H1I2|) ') that a logarithmic singularity at 
Lo = ve persists no matter what strength of impurity 
scattering is present. Thus, somewhat counterintuitively, 
impurity scattering does not smear out the singularity at 
uj = ve- As will be evident below, the singularity at 
LD = ve results in a collapse of the renormalized EDOS 
at that frequency. The phonon contribution to the imag- 
inary part of the self-energy is zero until lj = ve, and 
the impurity contribution is also zero aX lo ^ ve, due to 
the collapse of the electron band at this frequency. This 
combination ensures that the overall imaginary part of 
the self energy drops to zero at this frequency, as Fig. 3b 
shows. However, at w = 2ve, the result is in general 
non-zero: while the electron-phonon contribution is zero, 
since it depends on the EDOS with frequency shifted by 
ve, the impurity contribution is not zero since the renor- 
malized EDOS at that frequency is non-zero. Therefore 
the total function is non-zero for lo = 2ve- 

As Cappelluti and Pietronero pointed out already in 
Ref. HJ), the slope of the real part of the self energy at 
the Fermi energy {lo — 0) changes systematically with 
impurity scattering. This makes inference of the electron- 
phonon coupling strength from the self energy a very dif- 
ficult problem. Moreover, the "kink" observed in pho- 
toemission spectroscopy will be affected. However, 
for an Einstein spectrum, the "kink" actually becomes 
more pronounced with increasing impurity scattering, as 
shown in Fig. 4. This is because the logarithmic singu- 
larity remains in the real part of the self energy, but im- 
purity scattering makes this singularity sharper, as Fig. 4 
illustrates. As shown below, with a more realistic elec- 
tron phonon coupling spectrum the conclusion of Ref. 4 
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FIG. 5: The spectral function at the Fermi level (a) and the 
renormalized density of states for the parameters of Fig. 3, 
for various degrees of impurity scattering. Note that in (b) 
the "collapse" of the band evident at uj = ve remains for all 
strengths of impurity scattering. Spectral weight is pushed 
to higher energies as the impurity scattering increases. In 
(a) when no impurity scattering is present (solid curve) the 
spectral function is given by a delta function at the origin 
(indicated by the vertical line with the arrow), followed by an 
incoherent piece that starts aX lu = i^e = meV. 



applies, namely that the "kink" becomes smeared with 
increasing electron-impurity scattering. 

In Fig. 5a and 5b we show the spectral function at the 
Fermi energy and the renormalized EDOS, respectively, 
as a function of frequency. Note the significant reduction 
near a; = with increased impurity scattering in both 
properties, with the compensating weight gain occurring 
at high frequencies. 
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FIG. 6: Real (a) and imaginary (b) parts of the electron 
self energy, calculated using an interaction with phonons de- 
scribed by a Lorentzian spectrum centered at = 10 meV 
with A = 1. The Lorentzian describing the phonons has 
broadening parameter, 5 = 3 meV. Note that overall trends 
are the same as for the Einstein phonon (Fig. 3) except ev- 
erything is smooth. In the real part the degree of the dip 
(near u fa ve = W meV) is strongly altered by the impurity 
scattering (in the Einstein phonon case it remains singular). 



(B) LORENTZIAN PHONON SPECTRUM 

Many of these results contain anomalously sharp fea- 
tures because of the use of an Einstein model for the 

phonon spectral function. It was already shown in Ref. 

that these features are considerably smoothened when a 
more realistic spectrum is used. Here we use the broad- 
ened Lorentzian spectrum of Ref. with (5 = 3 meV. 
Results are shown in Fig. 6 for (a) the real and (b) imag- 
inary parts of the self energy. Clearly the singularity in 
the real part is absent even without impurity scattering; 
in contrast to the Einstein mode case, however, impurity 
scattering alters the nature of the dip near ve consider- 
ably. The rest of the frequency dependence for both real 
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FIG. 7: The dispersion corresponding to the parameters of 
Fig. 6. The severity of the "kink" at = 10 meV is signifi- 
cantly moderated by impurity scattering. 



and imaginary parts is also affected considerably, as was 
the case with the Einstein spectrum. In particular, the 
high frequency portion of the real part remains enhanced 
due to impurity scattering, and, as remarked in Ref. 0, 
this will affect the interpretation of the high frequency 
extrapolation of the dispersion relation as measured in 
photoemission [loj . Of special note here, however, is 
the definite reduction of the magnitude of the slope of 
the real part of the self energy near w = 0. The way 
this would appear for the dispersion relation is shown 
in Fig. 7. Note the ambiguity which impurity scattering 
causes for the extraction of an electron-phonon coupling 
strength from such a measurement 4] . 

Examining the entire spectral function and/or EDOS 
may be a more promising way to proceed. Of course some 
experimental challenges concerning background and res- 
olution would have to be overcome in the meantime. In 
Fig. 8 wc show the spectral function for (a) Cfc = and 
(b) efc = — 5 meV. Clearly once impurity scattering is 
included the separation of the quasiparticle contribution 
(normally a delta-function without impurities) from the 
incoherent part (normally with weight A/ (1-1- A) when the 
Fermi energy is infinite) will become very difficult. This 
is especially true for a broad electron-phonon spectrum, 
particularly one that goes to zero at zero frequency (the 
one used here is cut off before zero frequency is achieved, 
at w = 0.05meV). 

In Fig. 9, we show the renormalized EDOS, for vari- 
ous impurity scattering strengths. Even if the underlying 
electron band structure were known with absolute cer- 
tainty it would be very difficult to disentangle electron- 
phonon effects from electron-impurity effects. 

All of the results presented so far utilized a Lorentzian- 
shaped bare EDOS to represent the "bare" electron spec- 
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FIG. 9: The renormalized ED OS as a function of frequency for 
the parameters of Fig. 6. The bare EDOS is also shown. Note 
that the structure near the phonon frequencies is smoothed 
considerably by the impurity scattering. The latter also ac- 
counts for a considerable shift of spectral weight to higher 
frequency. 



V. SUMMARY 



-50 -40 -30 -20 -10 10 20 30 40 50 
(o(meV) 

FIG. 8: The spectral function at (a) the Fermi level and (b) -5 
meV for the same parameters as in Fig. 6. Note that a delta 
function occurs only at the Fermi level (indicated by the arrow 
in (a)) when no impurity scattering is used. Increasing the 
impurity scattering strength immediately broadens the delta 
function into a Lorentzian-like piece; the incoherent part is 
already broad, and changes only quantitatively. 



trum. The results in Ref. Q) were for a constant bare 
EDOS; these can be compared with one another to see 
what the effect of the bare EDOS is. Here we show a 
specific example in Fig. 10, for A = 1 and l/r = 10 meV. 
The low energy part tends to be very similar, while the 
higher energy features experience quantitative changes. 
This implies that low energy data will be somewhat in- 
sensitive to the bare band structure, whereas high energy 
data will depend quantitatively on the specific form of the 
band structure. These effects will make the analysis of 
ARPES data much harder as different effects are highly 
entangled. In addition, the high energy part of the vari- 
ous single particle properties discussed here is very much 
dependent on the underlying EDOS. 



In summary, we have examined the effect of both 
electron-phonon and single impurity scattering on the 
electron self energy and Green function under special 
circumstances. These include the presence of a finite 
bandwidth and phonons which are not infinitely sharp 
in energy. The broadened nature of the phonon spec- 
trum can be due to intrinsic effects on a single phonon 
(anharmonicity, interactions, etc.) that give rise to a 
broadened lineshape; on the other hand they may arise 
from simply accounting for the dispersion displayed by 
phonons in real materials. We found that the "smear- 
ing" effect expected from impurity scattering is much 
more pronounced when such a realistic phonon spectrum 
is used. Otherwise many single electron properties re- 
tain sharp features in frequency space when an Einstein 
spectrum is used. Unfortunately, this combination of fi- 
nite bandwidth (usually with unknown band structure) 
and uncertain degree of impurity scattering hinders the 
ability to infer details concerning the electron-phonon in- 
teraction. 
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FIG. 10: A comparison of the (a) real, and (b) imaginary part 
of the self energy, and (c) the dispersion, for a Lorentzian 
(solid curve) vs. constant (dashed curve) bare EDOS. Mostly 
high frequency parts change quantitatively. This means that 
a full inversion of such curves to determine the underlying 
microscopic interactions would in practice be very difficult. 



